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With the experimental data from STAR on the centrality dependence of transverse momentum 
Pt spectra of pions and protons in Au+Au coUisions at y^sjvjv ~ 200 GeV, we investigate the 
scahng properties of transverse energy Et distributions at different centralities. In the framework 
of cluster formation and decay mechanism for particle production, the universal transverse energy 
distributions for pion and proton can be described separately but not simultaneously. 
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I. INTRODUCTION 

Recently, the scaling behaviors of the transverse mo- 
mentum distributions for particles produced in high en- 
ergy collisions attract more and more attentions. The 
search on the scaling behaviors of particle spectrum is 
significantly important for understanding the evolution 
of high energy heavy ion collisions and particle produc- 
tion mechanisms, because the produced particle distribu- 
tion carries information about the dynamics of the sys- 
tem and is one of the most important observables in high 
energy collisions. In [H, with the experimental data 
from STAR, PHENIX and BRAHMS, we found there is 
a scaling behavior of transverse momentum px spectrum 
for pions. The scaling function is independent of the 
colliding system, the colliding energy, the centrality and 
(pseudo)rapidity. For protons and anti-protons, there 
also exists a similar scaling behavior which is indepen- 
dent of the centrality and rapidity in Au+Au collisions 
at Y^ijvw — 200 GeV For strange particles, such as 
K, A and S, the scaling behaviors for their distributions 
are under consideration. 

In high energy heavy ion collisions, a lot of the initial 
kinetic energy of the bombarding nuclei is converted into 
those of the produced particles, including longitudinal 
and transverse ones. Different from momentum, kinetic 
energy is a scalar and is directly associated with the tem- 
perature of the hot medium created in the collisions. For 
different species of particles, the same momentum corre- 
sponds to different kinetic energy because of mass effect. 
As well-known, the kinetic energy, rather than momen- 
tum, of particles in a thermalized system satisfies the 
Boltzmann distribution. Thus the distribution of kinetic 
energy of particles produced in ultra-relativistic heavy 
ion collisions is more effective in revealing the thermal 
properties of the system. With above consideration in 
mind, one can ask whether the transverse kinetic energy 
Et distribution has similar scaling behavior. This is our 
motivation of investigating the scaling properties of dis- 
tributions of transverse kinetic energy Et of particles. 
We hope the mass effect can be suppressed in the new 
scaling function. 



In this paper the scaling Et distribution of protons in 
the mid-rapidity region is studied and compared with the 
scaling pt distribution. It is very essential to ask why the 
scaling behaviors exist for different particles and what is 
the potential universal dynamics. In we found that the 
string overlap mechanism can not explain data for pion 
and proton simultaneously. In this paper we will consider 
another mechanism based on the parton percolation the- 
ory [1, @, 0], and hope to get more information about 
particle production mechanism in nuclear collisions. 

This paper is organized as follows: In Sec. II, we dis- 
cuss relation between distributions for Et and pT, and 
the scaling property of Et distribution of protons pro- 
duced in Au+Au collisions at ^/snn — 200 GeV. Uni- 
versal scaling functions for pion and proton are given. In 
Sec. Ill we discuss the universal scaling functions in the 
framework of cluster decay. The conclusion is drawn in 
Sec. IV. 

II. THE SCALING BEHAVIOR OF TRANSVERSE 
ENERGY DISTRIBUTION 

From the definition of transverse kinetic energy Et = 
rriT—mo = \/ Pt + fno—mo, with toq mass of the particle, 
the invariant Et distribution is 

rno^^fl^^ (1) 
ETdETdy Et pTdpTdy 

From the definition of Et , one can see that when mo — > 
or Pt ^ oo, Et is approximately equal to pT, and 

(2) 

ETdETdy pTdpTdy 

So it is feasible to only consider the scaling behavior of 
Pt distribution in the range pT > 0.5 GeV/c for pions, 
whose mass is only 0.139 GeV. Actually, the results ob- 
tained in P, 0, are excellent. But if mg is large, as for 
kaons, protons and anti-protons, etc, we must consider 
the mass effect and investigate the scaling behaviors of 
Et distributions in the same range of pt- 

In [il we have summarized the method for searching 
the scaling behavior of the particle's pt spectrum. Now 
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we can use the similar method to investigate the scal- 
ing behavior of the transverse energy distribution of final 
state particles, with transverse kinetic energy Et instead 
of transverse momentum px- First, define a scaled vari- 
able 



z = Et/K , 
and the scaled spectrum 

^ z) — A 

^ ' 2'KETdETd-q 



(3) 



(4) 



Et=Kz 



with A and K free parameters chosen to fit the scaled 
distributions to the same curve. Of course values of A 
and K are different for distributions at different centrali- 
ties, given that for most central collisions they are set to 
be 1. Then ^{z) can be regarded as a parameterization 
of the Et distribution in most central Au-|-Au collisions. 
To obtain a universal scaling function for all centralities, 
independent of the arbitrary in choosing values of A and 
K for the most central collisions, we introduce another 
scaled variable 



u = z/{z) = Et/{Et) 
and the normalized scaling function 



^'(u) = (z)2$((z)u)/ / ^{z)zdz 
Jo 

Here (z) is defined as 



z^{z)zdz/ / ^{z)zdz , 



(5) 



(6) 



(7) 



and one can easily check that J duu^j(u) — 1 and (u) — 
J duv?"^{u) = 1. With above steps, one can investigate 
whether there exists a scaling behavior of transverse en- 
ergy distribution for protons produced at mid-rapidity in 
Au-|-Au collisions at ^Jsnn = 200 GeV, which has a wide 
rang pT coverage Q . As shown in Fig. [U all data points 
for different centralities can be shifted to the same curve 
within errorbar, and the corresponding values of param- 
eters A and K shown in TABLE I. To parameterize the 
curve, we define v = ln(l + z), and the scaling function 
is 

$(z) = 18.604 exp(-0.6765u - 6.48^^ -I- 1.477w^) . (8) 

After normalization, 

= 1.4846 exp(2.394u - 5.79v^ + 0.977v^) , (9) 

with redefined v — ln(l + u). 

In order to see the agreement between the scaling dis- 
tribution and the data in linear scale, a ratio B can be 
defined as 

B=experimental data/fitted results. 
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FIG. 1: Scaling behavior of the Et spectra for protons pro- 
duced at mid-rapidity in Au-|-Au collisions at RHIC. Feed- 
down corrections are considered in the data. The solid curve 
is from Eq. (|9]). The data are taken from [i| after calculation 
with Eq. HJ. 
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0-12% 


0.4192 


0.0038 


0.5786 


0.1440 


10-20% 


0.4293 


0.0059 


0.5769 


0.1788 


20-40% 


0.4428 


0.0114 


0.5710 


0.3022 


40-60% 


0.4535 


0.0307 


0.5635 


0.6732 


60-80% 


0.4571 


0.1025 


0.5437 


2.0985 


40-80% 


0.4536 


0.0466 







TABLE I: Parameters for coalescing all data points to the 
same curves in Figs. [T]for p and |4] for tt^. 



First let us look at the agreement of the fit to the pt 
distribution. The fitted normalized pt scaling function 
of proton is obtained from [3], which is given as follows. 



0.064 exp(13.6'i; - 16.67w^ + 3.6v^) 



(10) 



with V = ln(l -|- u) and u = pt/{pt)- The result of B 
for the deviation of the scaling pT function from data in 
most central collisions is shown in Fig. [2l They agree 
with each other within experimental error. Then one 
can investigate B for the Et case. The result is shown in 
Fig. [3] As can be seen in Figs. [2]and[3l in the middle re- 
gion, both fits are good. But at the low and high region, 
the agreement of the scaling Et distribution is better. 
This indicates that the Et scaling can describe the ex- 
perimental data better. For data at other centralities, 
similar conclusion can be made. 
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FIG. 2; The ratio B of the experimental data to the scaUng 
Pt fitted results Eq. (|f Op . here u = Pt/{pt)- Errorbars shown 
are calculated from those from experiment with the definition 
of B. 
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FIG. 3; The ratio B of the experimental data to the scaling 
Et fitted results Eq. Q, here u = Et/{Et). 



III. SCALING TRANSVERSE KINETIC ENERGY 
DISTRIBUTION AND CLUSTER DECAY 

In and last section, we have found the scaling 

laws of transverse momentum and kinetic energy spectra 
for pions, protons and anti-protons . With those results, 
we should try to trace out the particle production mecha- 
nism. In [J] we have shown that string picture for particle 
production can not describe the scaling properties of px 
distributions for both pions and protons simultaneously, 
because from the pion and proton spectra one can get 
opposite changes of string overlap degree from central to 
peripheral collisions. 

In the framework of parton percolation [1, 0, 0| , mas- 
sive color-neutral clusters are formed in high energy colli- 
sions from produced partons. In this picture all observed 
final state particles are decay products of those clusters 



with different sizes. Now we discuss whether such a pic- 
ture can describe the transverse energy distribution. One 
can use the inverse of transverse kinetic energy squared, 
X, to describe the cluster size with distribution denoted 
by W{x). Here the transverse kinetic energy is used in 
place of transverse momentum in 0, 0, 0] ■ The normal- 
ized W(x) has been supposed to be a Gamma distribu- 
tion 



W{x) = :^{jx)'' iexp(-7x) 



(11) 



with 7 and k two parameters. To get the spectra of final 
state particles one needs to know the cluster fragmenta- 
tion function f{x,ET) for each species of particles. Here 
f{x, Et) gives the probability of producing a hadron with 
transverse kinetic energy Et from a cluster of size x. 
We have no first-principle as an instruction for the func- 
tional form of fragmentation function. We expect that 
the fragmentation function is a function of the fraction 
z = Et\/x of transverse kinetic energy of the produced 
particle relative to that of a cluster with size x. This 
fraction is a close analogue of that for usual fragmenta- 
tion functions from hard partons to hadrons. For this 
reason, one can assume that the fragmentation function 
for cluster decay takes the same functional form as the 
usual ones, as used in P] and references therein. 



f{z) = Dz%\-z)\l^cz'') 



(12) 



with D, o, 6, c and d five parameters. 

With the cluster size distribution and fragmentation 
function, the Et distribution of a species of final state 
particles can be expressed as 



dN 



c 



1/E- 



dxf{ETVx) X W{x) 



(13) 



In last equation C is the normalization constant for the 
total number of clusters formed before hadronization. In 
real fitting parameters C and D always appear as a prod- 
uct. So one can absorb D into C. The upper limit 
in the integration is given from the consideration that 
hadron's fraction of transverse kinetic energy cannot be 
larger than 1 . The physics behind above formula is as fol- 
lows. The clusters are formed before their decay, so that 
W{x) is the same for all species of final state hadrons. 
On the other hand, f{z) is different for different hadrons, 
because different hadrons come from different fragmen- 
tation channels of clusters, but f{z) should have no con- 
nection with the cluster size x, thus is independent of the 
colliding system, centrality and rapidity, namely univer- 
sal for all collision processes. This is also a parallel to the 
usual parton fragmentation functions which are indepen- 
dent of the processes for the hard parton production. It is 
easy to show that last equation guarantees scaling of the 
particle distribution. Under the transformation x —^ Ax, 
7 ^ 7/A and Et Et/V^, both the two functions 
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W{x) and f{z) are all invariant, so the Et distribution 
is also invariant within a normalization constant. This 
invariance means that a change of the mean transverse 
kinetic energy can be equilibrated by a change of the 
value of parameter A for the cluster distribution. Thus it 
shares the essence with the well-known renormalization 
group transformation. This invariance is the scaling we 
are looking for. 
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FIG. 4: Normalized scaling distribution for pions produced at 
mid-rapidity in Au-|-Au collisions at RHIC with scaling vari- 
able u = Pt/ (pt) ■ Feed-down corrections also are considered 
in the data. The solid curve is from Eq . (I13|l . and the dash 
curve is from Eq . (|14p . The data is from ^8j. 

Now we show how to get the Et distributions of pions 
and protons from the decay of clusters. Because pion is 
very light, Et — Pt in almost the whole range of observed 
transverse momentum, one only need consider the scaling 
Pt distribution. We use a recent set of data in [§] for 
pion. Considering pT range in the data used now is much 
wider than that used in [l| , we fit the new data. All data 
points can be put to the same curve by suitable A and 
K, whose values are tabulated in TABLE I. The new 
parameterization for the normalized scaling function is 
now 



(l + 13.6e-^-^"), (14) 



= 0.288 1 



7.868 



with u — pt/{pt)- As shown in Fig [51 the agreement be- 
tween the scaling function and data is great. To describe 
the scaling function for pion with the cluster mechanism, 
one can work with Eq. (fT3|) and fit the obtained scaling 
function with 7 parameters, C, 7 and k for cluster distri- 
butions, a,b,c and d for the fragmentation function from 
a cluster to pion. The fitted parameters are tabulated 
in TABLE II. For comparison, the fitted curve from the 
cluster decay is also plotted in Fig. IH Two curves for 
the scaling function and the fitted result cannot be distin- 
guished in the plot in log scale. For protons, one should 
consider the scaling Et distribution, as discussed in last 



section. To fit proton's scaling function, one cannot treat 
all the three parameters for the cluster distribution as 
free parameters. As will be addressed below, parameter 
k must be the same in fitting the scaling functions for 
pion and proton. The fitted parameters for proton case 
are also given in TABLE II, and the fitted curve is shown 
in Fig. [5] The normalized scaling Et function Eq .([9]) is 
shown in this figure too. As can be seen from Fig. [H the 
fit is also beautiful in a wide range of Et- The inset in 
Fig. [5] gives the ratio of the scaling function for proton 
to the fitted distribution from cluster decay. It should be 
noted that the mean transverse momenta are different for 
pion and proton in the collisions with given centrality. In 
obtaining the scaling functions we try to scale the mo- 
mentum (or transverse kinetic energy) to a variable with 
unit mean. For this purpose, the scaling factors for pions 
and protons must be different. The scale in pT or Et can 
be translated into that of x and finally to 7, as discussed 
in the paragraph following Eq. p3|l . So the values of 7 
for pions and protons obtained in our fitting are differ- 
ent. This difference in scaling factor also causes different 
values of C, as shown in TABLE II. Another reason for 
the difference of values of C for pion and proton is that C 
contains in fact parameter 13 as a factor which should be 
much smaller for proton than for pion because of the low 
yield of proton relative to pion. Besides, the distribution 
of clusters W{x) is the same for all species of final state 
particles, so one can figure out k is the same for pions 
and protons, because the scale change in pT or Et does 
not affect k. Because the parameter a for pion is about 
— 1 while that for proton is only about 0.2, the density 
of pion with low transverse kinetic energy (or momen- 
tum) is much higher than for proton, in agreement with 
experimental fact that p/n as pt 0. Now we 
can say the mechanism of cluster formation and decay 
can describe the universal scaling behaviors of transverse 
energy distributions for either pion or proton. 

But can the cluster mechanism describe all those final 
state particle spectra at the same time, as a true particle 
production mechanism should be able to? The answer 
lies in the trends of change of the parameter K for pion 
and proton from central to peripheral collisions. K in- 
creases for pions but decreases for protons. The same 
results are also obtained in earlier works P, 0, d, 01 for 
the Pt distributions. Such trends are closely associated 
with the fact that the pion spectrum is suppressed from 
intermediate to high pT in central Au-t- Au collisions while 
no suppression was observed for the proton spectrum. In 
fact, from Eq. ^ one can see that the values of K used 
in shifting data points to the normalized scaling distri- 
bution are the mean values of transverse kinetic energy 
(Et)- With the cluster mechanism, on the other hand, 
(Et) can be calculated from Eq. (flSll as 



K = {Et) = 



dxW(x)/x^/^ /q dzz^fjz) 

j;;°dxw(x)/x j^'dzzfiz) 



(15) 
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tons in Au+Au collisions at 200 GeV. We showed that 
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FIG. 5: Comparison between the normalized scaling distribu- 
tion for protons produced at mid-rapidity in Au-|-Au collisions 
at RHIC and the fitted result from the cluster decay mecha- 
nism. 



As stated above, cluster distribution W{x) is the same 
for different species of particles, while cluster fragmenta- 
tion functions f{z) are the same for a species of particle 
produced, independent of colliding centrality. If we take 
a ratio between values of K for different centralities, the 
ratio must be the same for all particles, since the f{z) 
terms cancelled in the ratio. This demands that values 
of K for different species of particles must change in the 
same way with centrality. This demand cannot be sat- 
isfied in the cluster mechanism. This fact indicates that 
the cluster formation and decay mechanism is not a uni- 
versal one for the production of final state pions, protons, 
and other particles in high energy collisions. As shown 
in 0, however, one can describe spectra of both pion 
and (anti)proton with different values of the parameter k 
in the cluster distribution. Different values of k implies 
that pions and (anti)protons originate from clusters with 
different distributions. 

Finally a brief discussion on the distributions at LHC 
energies can be addressed. In Pb+Pb collisions at LHC 
energies (Et) should be much larger than that at RHIC 
energies. From the universal distribution discussed in 
this paper, the value of A in the cluster distribution must 
be much smaller, which means that the clusters formed in 
those collisions at LHC energies must have mean trans- 
verse kinetic energy much larger than that at RHIC en- 
ergies. 



IV. CONCLUSION 

In this paper we discussed the scaling properties of the 
transverse kinetic energy distributions for pions and pro- 
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12.7766 


0.3857 


7 


203.7292 


318.4096 
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4.3164 


4.3164 


a 


-1.0842 


0.1890 


b 


10.0671 


25.9110 


c 


201.0009 


142.9827 


d 


0.0016 


0.5406 



TABLE II: Values of parameters of the universal transverse 
energy distributions Eq. (|13|l for pions and protons. 



better than that for the pr case. Although the cluster 
formation and decay mechanism can describe spectra of 
either pion or proton at different centralities, it cannot 
give consistent description, with the same parameters in 
the cluster distribution, for both species of particles at 
the same time. Thus it can be excluded as a consistent 
particle production mechanism. Other particle produc- 
tion mechanism must be in effect for the scaling behavior. 
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